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This paper describes the modeling of quantum dots lasers with the aim of assessing the conditions 
for stable cw dual-mode operation when the mode separation lies in the THz range. Several pos- 
sible models suited for InAs quantum dots in InP barriers are analytically evaluated, in particular 
quantum dots electrically coupled through a direct exchange of excitation by the wetting layer or 
quantum dots optically coupled through the homogeneous broadening of their optical gain. A stable 
dual-mode regime is shown possible in all cases when quantum dots are used as active layer whereas 
a gain medium of quantum well or bulk type inevitably leads to bistable behavior. The choice of a 
quantum dots gain medium perfectly matched the production of dual-mode lasers devoted to THz 
generation by photomixing. 
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I. INTRODUCTION 

CW THz and millimeter wave generation using the 
beating frequency of a dual-mode laser is a topic that has 
attracted many efforts during recent years because there 
is currently a lack of versatile and easy to use sources 
at these frequencies^. To fill this gap, THz photomixing 
has many strengths, including the use of photonic tech- 
nologies as the manufacture of laser diodes and complex 
optical systems or the ability to carry the beat on opti- 
cal fiber. In all cases the conversion of the optical signal 
is provided by a photomixer which transforms the fast 
variations of the optical intensity in the electromagnetic 
radiation. Since a few years these components exist at 
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once in GaAs technology for wavelengths 
InGaAs technology around s» 1.55 ^m 3 -. 

Regarding the optical beat source, several alternatives 
were considered. The easiest way is to combine two inde- 
pendent lasers that are stabilized in frequency and ampli- 
tude. This solution is now commercial but it is expensive 
because of the use of extended cavity lasers or high power 
DFB requiring sophisticated control electronics^. In ad- 
dition, the use of two independent lasers implies that 
the drift and noise are added in the beat noise which 
directly blames the linewidth of the THz signal gener- 
ated. The implementation of this solution in the case of 
high-resolution spectroscopy requires a mctrological sta- 
bilization of lasers that remains difficult in practice^. 

The use of a single laser seems more simple and should 
allow an improvement of compactness. Two solutions are 
possible. In the first the laser operates in pulsed mode 
either by gain-switching or by active or passive mode 



locking. The THz beating is then obtained by filter- 
ing the high-frequency harmonics of the pulse repetition. 
Because the consistency of the optical phase is main- 
tained in these modes, extremely fine laser linewidths 
can be obtained^. Besides tunability weakness in the 
case of mode-locking, the major drawback of this tech- 
nique lies in the sophisticated filtering process which is 
itself a source of high attenuation and that needs to be 
offset by several optical amplifiers. This undermines the 
simplicity and cost of final assembly. Alternatively the 
use of a suitable structure in the active layer allows for 
the dual-mode filtering task within the passively mode- 
locked laser itself. Millimeter frequency generation has 
been demonstrated that wayS, however, the complexity 
is transferred to the laser fabrication that must include 
a specific longitudinal structure along the optical waveg- 
uide. 

The second solution to generate THz with a sin- 
gle component is to use a laser operating in the dual- 
mode regime. Again several techniques have been pro- 
posed. First with the advent of manufacturing tech- 
nologies, semiconductor lasers can now include a com- 
plex structure with multiple DFB sections operating 
independently 9 '—. The fact that the two modes inject 
each other can lead to instabilities and may seriously de- 
crease the tuning range. An elegant solution has been 
demonstrated with a diode-pumped solid-state laser us- 
ing either a separation of the two polarization eigenstates 
or a spatial separation of the two modes using an intra- 
cavity birefringent crystal. These latter structures have 
shown excellent performance both in terms of noise and 
beating tunabilit y 12 ' 13 . 
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To simplify, we may want to use a single Fabry- 
Perot edge-emitting laser. Under certain conditions, it is 
then possible to obtain structures exhibiting a dual-mode 
beating at the roundtrip time of the optical cavit y 14 ' 15 . 
This was used once to build an very compact imaging sys- 
tem at millimeter waves 1 ^. Similarly we may want to use 
a vertical-cavity semiconductor-laser which by construc- 
tion has much less longitudinal modes available in the 
gain curve. One must then stack two cavities within the 
structur o 17 ' 18 to obtain the required two separate longitu- 
dinal modes spaced in frequency by an amount dependent 
of the optical coupling. This concept has been recently 
taken over by some of us with a single vertical-cavity 
Bragg-mirror resonator and a photonic-crystal support- 
ing a slow optical Bloch mode^ 9 -. These structures have 
many advantages for planar applications. Encapsulating 
the two cavities within the same small component ensures 
that any optical-cavity frequency-drift, for example due 
to temperature variations, will apply similarly to both 
modes. As a consequence the beating frequency will be 
insensitive to such drifts at first order. A built-in dual- 
mode laser is also an integrated component extremely 
compact and easily integrated into a system for which the 
stability of the beat will be easily obtained by monitoring 
the overall operating parameters of a single component 
(temperatures, currents, ...). To enjoy these benefits, 
however, it is necessary that the modal competition does 
not destroy the desired dual-mode operation. 

Modal competition is at the heart of the stable dual- 
mode operation of a single semiconductor laser. Since 
the seminal work of Lamb 2 ^, it is well known that even 
if the laser cavity allows two simultaneous modes with 
similar quality factors and gains, the gain medium non- 
linearity sometimes forbids their simultaneous emission. 
A stable dual-mode operation is only allowed when the 
coupling factor between the modes is not strong enough. 
Otherwise the laser operates in a bistable regime where 
one mode dominates. Although unwanted, this case un- 
fortunately happens with gain media such as bulk semi- 
conductor or quantum wells (QWs). The physical ori- 
gin is the short conduction-band intraband relaxation- 
time that couple two adjacent modes sharing the same 
carrier population. Following Agrawai 2 ^- and taking an 
intraband relaxation time below 100 fs which is consis- 
tent with InGaAsP active layers^ 2 - or InGaAs QW— , it 
is not possible to have a stable dual-mode operation if 
the two modes are separated by less than 1.6 THz. This 
result was also obtained from rate-equation o 24 ' 25 . Such 
an estimation moreover neglects the spatial hole burn- 
ing which is also involved in the coupling of longitudinal 
modes 2 ^. This leaves no chance for a semiconductor laser 
emitting at 1.55 /im and incorporating a bulk or QW gain 
medium to operate dual-mode with a beating of « 1 THz. 
The dynamics of multimode VCSEL 2 - or edge-emitting 
lase r 28 ' 29 , and even multisection DFB lasers 30 claims this. 

Our dual-mode laser design involving a slow optical 
Bloch mode in a vertical cavity structure is supposed 
to operate dual-mode^. It is thus mandatory to incor- 



porate a gain medium that will not preclude the laser to 
operate in the desired regime because of a too strong cou- 
pling between modes. As shown by the above analysis, 
this coupling is mainly due to the sharing of the same 
carrier populations in the conduction band by the two 
adjacent modes. The proposed solution to circumvent 
this strong coupling is to use InAs quantum dots (QDs). 
The purpose of this work is then to verify that either 
optical or electrical coupling between InAs dot popula- 
tions will remain at a sufficiently low level to maintain 
the dual-mode operation. Dots are optically coupled by 
the homogenous broadening of the optical gain that is 
estimated a few meV at 300 K^ir— . Dots are also elec- 
trically coupled by their common interaction with the 
wetting layer that may allow a transfer of excitation from 
one dot to another. 

The rest of the paper is organized as follows. Section 
ITT1 reminds the analytical modeling of dual-mode lasers 
with a single gain medium. Section ITTT1 details the case of 
dual-mode lasers including QDs by successively assum- 
ing first that dots are not coupled, and then coupled by 
a chemical-like equilibrium law linking excited dot pop- 
ulations through an interaction with the wetting layer, 
and finally assuming that the QDs are optically coupled 
through their homogeneous emission width. The conclu- 
sion is given in section llVl 

II. DUAL-MODE SEMICONDUCTOR LASER 
WITH A SINGLE GAIN MEDIUM 

A steady dual-mode operation requires special condi- 
tions for the allowed modes of a laser. The analysis ini- 
tially given by Lamb 2 ^ is still the most relevant and eas- 
ier to understand, it is based on two coupled differential 
equations that describe the time evolution of the modes 
in the laser cavity. These modes are characterized by 
their respective intensities I\ and I 2 

^± = (a 1 -p 1 I 1 -e 12 I 2 )I 1 (la) 
at 

-f = (a 2 - foh ~ 62J1) h (lb) 
at 

where a* are the unsaturated optical gains and where fii 
and Oij are the self-saturation and cross-saturation gain 
coefficients. The stable solutions of these equations are 
obtained by canceling the derivatives 

= ai - (8 12 /p 2 )a 2 ps = a 2 - (6>2i//3i)ai , . 
1 (l-C)fr 2 (l-CJft {> 

with 

_ 012#21 

a dimensionless coupling factor introduced by Lamb who 
showed using a perturbation analysis^ that it controls 
the existence and stability of the dual-mode operation. If 
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C > 1, i.e. in the case of a strong coupling between the 
modes, only single mode is enabled, the system is bistable 
and the laser operates continuously only on either of the 
two modes. 

Starting from this conditions, it may be possible how- 
ever to obtain a stable dual-mode regime changing the 
gain medium itself. This has been shown in several ways: 

• by adding a saturable absorber in the cavity 9 ; 

• by physically decoupling the two gain areas (for 
example by writing a dual-period DFB within the 
cavit y 36 ' 37 , or a solid-state laser where the gain is 
physically separated into two different paths, for 
example by changing the polarization^). 



Neither of these solutions is satisfactory and feasible 
for our application where both modes share the same gain 
location^. We will further investigate the link between 
this model and the most common rate-equations model 
of semiconductor lasers. 

A semiconductor laser with a bulk or QW gain medium 
is generally described on the basis of population rate evo- 
lutions. Ahmed and Yamada^ gave a fairly complete 
version of it from which one has extracted the equations 
describing a dual-mode laser 

^ = k - — - AxNSi - A 2 NS 2 (3a) 
at r n 

^ = -—+A 1 NS 1 (l-s 1 S 1 ~$ 1 S 2 ) (3b) 
at r Pl 

— - = + A 2 NS 2 (1 - e 2 S 2 - 2 Si) (3c) 

where N is the carrier density beyond transparency, 
k = I/(eV) is the pumping parameter which depends 
directly on the injection current I and volume of the ac- 
tive region V, r„ is the carrier recombination time, the 
Ai are the differential gains and Si the corresponding 
photon densities, r Pi are the photon lifetime within the 
cavity, £j and $i the self-saturation and cross-saturation 
coefficients of optical gains. Indices i refer of course to 
the various modes allowed by the cold cavity here reduced 
to 2. Considering a laser running well above threshold, it 
is possible to neglect the non-radiative recombination of 
carriers N / r„ whose contribution can be seen as implic- 
itly counted by a slight modification of the n- values which 
reflects the pumping threshold. In Eqs. ([3]) the coupling 
terms between the two modes appear explicitly with 
As their physical origin is similar to the self-saturation 
gain terms £j, these must also be taken into account. The 
omission of the correction factor (1 — EiSi — $iSj) in the 
equation carriers however provides a great simplification 
of the calculations at the cost of a very weak approxima- 
tion. 

Like with the Lamb model, simple calculations allow 
to obtain the stationary solutions of Eqs. Among 
the solutions, only the one with simultaneous nonzero 
optical intensities Si and S 2 corresponds to a dual-mode 



laser, others accounting instead of bistable operations. 
In the general case these solutions are quite complicated 
but in the simplistic case where the two modes are close 
to the gain maximum with a spacing small enough that 
their parameters are very similar, then £j =£,'!?, = t?, 
t P4 = t p and Ai = A and the stationary values are 



A 



(4) 



gss gss 

1 2 ~~ 2u 

with u = |(e + i?)k + 

The stability of this solution is further demonstrated 
by an analysis of variations around the operating point. 
To do this we replace in (j3|) rates by their stationary 
values of ([J} plus a small variation, N = N ss + n, 
Si = Sf s + Si. Keeping only the first order perturbations, 
we obtain a matrix system of three coupled first-order 
differential-equations that one writes in matrix form. 
The dual-mode stationary solution is then stable if the 
variations diminish over time, which is true if and only 
if the real parts of the eigenvalues of the corresponding 
matrix are negative. 

In our particular case, the matrix is 



/Are 







Ak 

Ak 
2u 2 r, 




(5) 



whose characteristic polynomial is 



V(X) 



A 



U \T p 2 



(6) 



The eigenvalues of are the roots of V, that is § ($ — e) 
and the two roots of the polynomial of degree two on 
the right side. With our choice of parameters, all coef- 
ficients of this polynomial are positive and then the two 
roots always exhibit negative real parts. Consequently, 
the dual-mode stationary solution proposed is stable if 
and only if $ < e. 

The self-saturation and cross-saturation coefficients 
can be calculated using a quantum model of the gain 
medium and the result shows typically d ~ 4g/3 21 i 38 ' 39 . 
In such conditions, natural dual-mode behavior is not ex- 
pected from a semiconductor laser, unless strong disper- 
sive effects are included in the cavity to remove the condi- 
tion assumed here: Si — e, ??; = r Pi — t p and Ai = A. 
This was recently confirmed by numerical simulations of 
the dynamics of semiconductor lasers made with typ- 
ical parameters from GaAs or InGaAsP technologies. 
Bistable behavior or highly multimode with a total power 
usually spread over a large number of modes is always ob- 
tained, in perfect agreement with experiments 2 ^^. This 
result was also confirmed with a Maxwellian modeling of 
laser— which only shows possibilities of either a single- 
mode behavior (actually bistable), or highly multimode. 



4 



Both our analysis and these works indicate that dual- 
mode semiconductor lasers are of course totally unex- 
pected! 

A similar extensive stability analysis was proposed ir>22 
to account for the degeneracy of the two polarization 
modes of VCSELs. This question is mathematically anal- 
ogous to our problem invoking two longitudinal modes 
except that there is no theoretical estimations of the cou- 
pling factor in that case. A fairly comprehensive system 
behavior was reported, including the polarization switch- 
ing of a VCSEL based on the gain saturation in a manner 
consistent with what is observed experimentally. 

III. QUANTUM-DOTS DUAL-MODE LASERS 

QD lasers are significantly different from bulk or QW 
semiconductor lasers discussed above. Indeed QDs once 
created have little interaction with each other and sub- 
stantially the same ability to capture a free carrier from 
the barrier. Several things differentiate strongly this type 
of material gain with previous ones: 

1 . The total number of active QD in resonance with a 
given mode is determined by the construction, and 
even if this number is high, it remains much smaller 
than the possible number of states allowed in the 
conduction band of a bulk or QW gain medium. 

2. QDs are much less interconnected than the ex- 
cited states in the conduction band of a more usual 
semiconductor laser. In particular, they cannot di- 
rectly exchange carriers either within the QD pop- 
ulation of the same mode, or between QD popu- 
lations addressing different modes. In practice the 
exchange of carriers between two QDs must involve 
first thermionic emission to the barrier and capture 
of this carrier by another QD. This indirect process 
considerably weakened the coupling as compared 
to the case of bulk or QW lasers. Given the values 
??of the energy levels at stake in a system with InAs 
QDs in InP barrier, we neglect this coupling here 
although it could be taken into account in a subse- 
quent Monte Carlo numerical model of dual-mode 
operation of the laser. 

3. According to the QDs manufacturing method, a 
monoatomic InAs wetting layer appearing at the 
growth interface is often common to all dots. Ow- 
ing to this layer, the QDs are likely to have a kind of 
direct electrical coupling that can transfer the exci- 
tation of a dot family to another family addressing 
another mode. The coupling effectiveness obviously 
depends on the average distance between dots, and 
on the nature and thickness of this wetting layer 
which is intimately dependent of the growth tech- 
nology. 

4. The homogeneous optical linewidth of a single 
QD is estimated between 3 and 12 meV at room 



temperature^—. A direct optical coupling thus 
necessarily occurs between two QD populations ad- 
dressing optical modes separated energetically by 
ps 4.2 meV to obtain a beating at ~ 1 THz. 

We will evaluate these situations in sequence and esti- 
mate analytically how the various coupling effects inter- 
act with the stability of the dual-mode laser. 

A. Uncoupled quantum-dots 

Let us first assume an ideal dual-mode QD laser with 
two uncoupled dots families perfectly centered on the two 
optical modes allowed by the dual-mode cavity. Rate- 
equations thus write 

^£ = aiK _ El. - Ai (Ni(Si + 1) - PA) (7a) 

OX T n i 

^ = -—+A i (N i (S i + l)-P i S i ) (7b) 
at t Pi 

with i = 1, 2 and E>i = Ni + Pi the total number of QD 
addressing mode i, Ni the number of excited dots and 
Pi the number of unexcited dots, k the total pumping 
and < ai < 1 with ot\ + 012 = 1 two parameters that 
reflect the distribution of the pump between the two QD 
families. As above Ai terms, r ni and r Pi account for the 
modal gains, the carrier lifetime in the excited state and 
the photon lifetime within the cavity for each of the two 
modes. 

The reality of a QD laser is not exactly the one de- 
scribed here. In practice, the size dispersion of quantum 
dots produces an inhomogenously broadened gain cover- 
ing almost uniformly the two allowed modes with respect 
of our assumption of a THz-beating dual-mode laser. As 
a result many dots are excited but unable to feed a laser 
mode, thereby contributing to an increase of the laser 
temperature and a degradation of its quantum efficiency. 
Although neglected in Eqs. (|7|), this effect could be taken 
into account by an adequate reduction of the effective 
pumping term k for both QD populations. 

A careful inspection of Eqs. ([7]) clearly shows that there 
is no direct coupling between the two QD populations. 
The stationary solution obtained by canceling the deriva- 
tives and keeping only the dual-mode physical solution 
with positive populations is 




where i = 1,2 depending on the considered mode and 
n = -j- + BiT Pi + T ni and t[ = 2a l Kr Pi T„j. Of course the 
number of unexcited QDs directly derived, P* s — B>i — 
Nf s . Similarly the threshold is calculated analytically 
as the pumping level where a break occurs in the slope 



5 



of the photon stationary solution, t[ = r, — 2r ni 4£ In 
fact, Eqs. © describe a behavior of two completely inde- 
pendent lasers with quantum efficiencies and thresholds 
evidencing no modal competition. Assuming Sf s 3> 1, it 
is easy to verify with the method proposed by Lamb that 
this solution is stable and that it can be generalized to 
any number of n juxtaposed modes. 



B. Quantum-dots electrically coupled through the 
wetting layer 

To go beyond this first simple view of a QD laser 
where the decoupling between modes is total and there- 
fore where no modal competition will affect the two-mode 
regime, we now imagine that the two QD populations can 
interact with each other directly. Presumably an excited 
dot of the first family can transfer to an unexcited dot of 
the second family, the difference in photon energy being 
supplied or absorbed by the states involved in the wetting 
layer. This direct interaction is similar to a chemical-like 
equilibrium law between populations 



Ni + P 2 O N 2 + Pi 
In practice Eqs. (O becomes 



dNi 
~~dT 

dN 2 
dt 

dSi 
~dT 

dS 2 
dt 



Ni 

= axK + k 1 N 2 Pi - k 2 NxP 2 - 



PiSx 



a 2 K + k 2 NxP 2 - fciiV 2 Pi 



N-, 



A 2 (N 2 (S 2 + 1) - P 2 S 2 ) 
Si 

T, 



'm 
T p 2 



+ Ai (N 1 (S 1 + 1)-P 1 S 1 ) 
+ A 2 {N 2 {S 2 + 1) - P 2 S 2 ) 



(9) 

(10a) 

(10b) 

(10c) 
(lOd) 



where fc, are constants reflecting the exchange rate of 
the excited states between the two QD families. Clearly 
if fcj — QDs are not coupled anymore and one retrieve 
the previous situation. Instead, if ki constants are strong 
enough their corresponding rates in Eqs. (|10[) quickly be- 
come of paramount importance as compared to any other 
rate. This is because they are assigned to products of two 
carrier populations that correspond to the highest values 
involved in Eqs. (|10[) . Thus the proper operation of the 
laser is obtained when these two rates almost annihilate 
or equivalcntly 



k 2 NiP 2 a /ci7V 2 Pi 



h 

NtP 2 « CN 2 Px with C = — 

k 2 



mode in Eqs. (jT0|) . Dual- mode stationary solution for 
carriers is then simply given by 



2 



1 



with i = 1, 2 



(ii) 



But the dual-mode stationary solution for photons is 
much more complicated in the general case. Two spe- 
cial cases are however interesting to consider: 

(a) when the coupling induced by the wetting layer are 
identical, ki = k 



'Pi 



1 / B * T P t 



AiT n , 



+ kr Pl 



Ai T Pi 



B, 



AjTpj 



(12) 



(b) when the material and cavity parameters are iden- 
tical for the two modes, i.e. A, t = A, Bi = B, 



Sf = CtiKTp 



B 2 



2A 2 i 



(13) 



where — (1,2) or (2,1). If cases (a) and (b) occur 
simultaneously, the situation is similar to that of ^III Al 
because the /c^-dependent terms cancel in Eqs. (TT2"j) and 
(TT3")) . therefore the coupling by the wetting layer becomes 
inefficient. Out of this particular case, the coupling by 
the wetting layer induces a differentiation between the 
thresholds of the two laser modes obtained approximately 
here by setting Si = in (12]) et [fig]). The mode that had 
the lowest threshold before coupling feeds its population 
of excited QDs by direct transfer from the one who had 
the highest threshold, the result is an amplification of 
the threshold difference. In the extreme case where the 
constants ki are larges, the population of excited QDs 
which reaches the first its stimulated emission threshold 
is clamped and in turn saturates the population of excited 
QDs for the other mode to a level below its own threshold. 
Strong coupling between QD populations thus leads to 
the destruction of the two-mode regime. 

In the case where the coupling is sufficiently moder- 
ate so that the dual-mode laser regime is preserved, the 
stability analysis is conducted as in fJTX| We report the 
stationary values plus a small fluctuation in the rate 
equations ([TO]) to obtain the evolution matrix of these 
fluctuations 



which is the analogue of chemical equilibrium described 
by Eq. ([§]) where carrier populations are strongly coupled. 

For the sake of simplicity we now neglect the terms in 
this '+1' that account for spontaneous emission in each 







-ax Q 2 -<5i 

Ci ~ °2 — S 2 

Px 

o p 2 



(14) 
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with 

01 = 2A 1 ST 2 = 2A 2 Sf 

Ci = hN s 2 s + fc 2 P 2 ss C2 = feivf + fciPf 

5i = A x {Nf - Pf) 5 2 = A 2 (Nl s - P 2 SS ) 

ci = 0i + Ci + 1/ T m c 2 = ^2 + C2 + l/r» 2 

which are all positive terms since the population inversion 
above the threshold requires Nf s > Pf s . 

Fluctuations may return to equilibrium if and only if 
the eigenvalues A of have negative real part. These 
eigenvalues are the roots of the characteristic polynomial 

P(A) = A 4 + (0-1 + o- 2 )A 3 + (7T1 + 7T 2 + 0-10-2 - CiC2)A 2 
+ (7110-2 +7T 2 0-l)A + 7Ti7r 2 . (15) 

Because o^ > Q, all the coefficients of P(A) are posi- 
tive. According to the Routh-Hurwitz criterion, all roots 
will then have their real part negative if and only i^i 

((71 + <7 2 )(7ri + 7T 2 + CTlO-2 - ClC 2 ) 

/ , ^ ^ ( ai + cr 2) 2 7ri7r 2 

- (7TiCr 2 + 7T 2 Ci) > (16) 

7I"lC 2 + 7r20"i 

We can rewrite the difference of the two members of (fTBT) 



(Ci + 0- 2 ) (0"i 0- 2 - Cl C 2 ) + (""1 0-1 + 7T 2 0- 2 ) - 



(a 1 + 0-2) 2 7Tl7T2 

7Tl0 2 + 7T 2 0-1 

/ ^ \ 2 

= (0-1 + o- 2 )(oio- 2 - C1C2) + oio- 2 — — - — (17) 

~K\(J 2 + 7T 2 0-1 

and remark that since o"j > Ci this difference is always 
positive. We then conclude that the dual-mode station- 
ary solution obtained from Eqs. (fTOf is stable. It proves 
that the coupling of carriers by a chemical-like equilib- 
rium law between the two QD populations does not pro- 
hibit a dual-mode behavior in the sense of modal com- 
petition, even if it is likely to push the threshold up to 
an unacceptable value for one of the lasing mode if too 
strong. 



C. Quantum-dots optically coupled through the 
homogeneous linewidth 

We now consider the two QD families coupled only 
through the homogeneous gain width of each dot, the 
rate equations then become 

dN-i Ni 

^=a 1 K-^--A 1 (N 1 -P 1 )(S 1 + eS 2 ) (18a) 
at t„i 

dNn No 

_L = a2K _ _2 _ A2 ( N2 _ p 2 )( 52 + e5l) (18b) 
at r„9 



^ = -— + A ± (Ni - Pl)Si + eA 2 (N 2 - P 2 )Sx 
at t p , 

(18c) 

= -— + ^ 2 (iV 2 - P 2 )S 2 + eAxiNx - P 1 )S 2 



with of course the particle conservation of each QD family 
N~i + Pi = Bi and assuming above threshold operation for 
each mode which removes the spontaneous emission term 
'+1' in photon density equations. The variable e here ac- 
counts for the direct optical coupling between dots of the 
two modes. It depends on the value of the homogeneous 
gain width and thus ranges from « 10% to 60% for a 
1 THz frequency separation, depending on whether one 
considers a broadening of 3 meV or 10 meV. 

In the general case the expression of the stationary 
solutions is complex, but in the particular case of two 
populations having the same QD material parameters 
Tni = T n , T ni — t p , Ai — A, Bi = B, and identical 
1/2, the stationary solutions are 

1 



' ni 1 r. 

pumps a 



N ss = — 4 
2 

ess _ KT P 



2A{1 
Bt„ 



1 



(19a) 
(19b) 



2r„ 2A(1 + e)r„ 

These formulas are close to those obtained in mil AI for 
uncoupled QDs. In fact only the term 1 + e differs and 
makes the calculated thresholds apparently lower in this 
case as compared to uncoupled QDs. This distortion 
effect does not hold if we consider that the differential 
gain of a laser with high homogeneous linewidth will in- 
evitably be far lower than that of a laser with a narrower 
linewidth. No decrease in thresholds related to modal 
coupling should be expected as a simplistic reading of 
Eq. (|19b|) suggests. 

Like in previous section with carriers coupled by the 
wetting layer, the stability analysis of the dual-mode 
regime has been conducted in the general case. Again 
this stability is governed by the sign of the real parts of 
the eigenvalues of the evolution matrix 



= 








-<*1 


-eo"i 





-c 2 


-eS 2 


-s 2 


01 










6 


02 









(20) 



Again we have chosen to define as positive all the terms 
of this matrix 



0i = 2AiS? 

Ci = 2e^i^ s 

5x =A 1 (N 1 BS -P^) 

ci = 0i + £1 + 1/7-ni 



02 = 2A 2 S S 2 S 

C 2 = 2eA 2 ST 

02 = A 2 (7V 2 SS -P 2 SS ) 

C2 = 02 +£2 + l/r n2 



(18d) 



and the eigenvalues are the roots of the following charac- 
teristic polynomial (assuming 7Tj = <5i(e£i + 0i) for further 
simplification) 

P(A) = A 4 + (01 + a 2 )A 3 + (0-10-2 + Tn + tt 2 )A 2 

+ (TrKX, + 7T 2 ai)A + (1 - e 2 )<5i<S 2 (0i0 2 - (21) 

As £i£ 2 = e 2 0i0 2 , the constant term can be written 
more simply as c = (1 — e 2 ) 2 <5i<5 2 0i0 2 , and all the coeffi- 
cients of P(A) are then positive since < e < 1. Accord- 
ing to the Routh-Hurwitz criterion, P(A) roots exhibit 
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negative real parts if and only if the following amount r 
is positive^ 



r = (0'l+0'2)(7ri+7r2+0'l<72)-(7ri0'2+7r2<Tl) 

7Tl (T 2 + 7T2CT 1 

yielding after simplification 
r = (oi + CT2)cri(72 + (ttiO"i + 7r 2 (T2) — ° X ~*~ — — — (22) 

7Tl(T2 + 7T2(Tl 

where we see that r is an increasing function of e on [0, 1] 
as a sum of functions of such type. The required proof 
is thus obtained by just checking if the value at e = is 
positive. Noting that at e = 0, c = 7Ti7T2, allows to write 
as in JTIllEq. (H7J) 

, s (cti+ct 2 ) 2 c (7Ti-7r 2 ) 2 

{Hl<J\ + 7T2<72j = CT102 

7Ti<72 + 7T2<7i 7Ti (72 + 7T2 <J\ 

This proves that r > and concludes again to the sta- 
bility of the dual-mode regime in that case. It should be 
noted, however, that if e tends to 1, which corresponds 
to a gain width very large compared to the energy gap of 
the two modes, then two eigenvalues tend to zero making 
the behavior of the dual-mode laser extremely unstable. 
This is a limit case where the laser is very to close to 
a bistable behavior but this is not the scenario we have 
chosen for THz radiation by photomixing with an energy 



gap of 4.1 meV between modes and QD homogeneous 
linewidths between 3 and 10 meV. 



IV. CONCLUSION 

We studied different semiconductor laser types devoted 
to operate continuously in stable dual-mode emission. 
The proposed analysis is based on the rate-equations and 
Lamb's theory. After considering the case of laser with 
bulk or QW gain mediums for which it is demonstrated 
that the natural behavior is bistable and not dual-mode, 
we evaluated analytically what should be expected from 
QD semiconductor lasers. We have considered uncoupled 
as well as coupled QDs either electrically through a direct 
exchange of excitation by the wetting layer, or optically 
through the homogeneous broadening of the gain. In all 
cases we have shown analytically that a stable dual-mode 
emission is possible. 

If the objective of building a dual-mode semiconduc- 
tor laser producing a stable beating frequency in the 
THz range is highlighted, the selection of a QDs as gain 
medium is the most suitable as compared to bulk or QWs 
that inevitably lead to a bistable behavior. The incorpo- 
ration of QDs in the active membranes of future photonic 
crystal optoelectronic components for the 2.5D THz ra- 
diation by photomixing is on the way^ 2 -. 
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